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According to the scaling theory of localization, all quantum electronic states are localized in two- 
dimensional (2D) systems. One consequence of the theory is that there is no quantum percolation 
transition in 2D. However, in a real system at a finite temperature, electron phase coherent length 
is finite and the system is between quantum and classical. We find, in such a 2D system, a metal- 
insulator transition (MIT) caused by a novel type of percolation, meta-percolation. The relevance 
to recently observed 2D MIT is also discussed. 

PACS numbers: 71.30.-|-h, 73.40.Hm 



There have been long lasting interests to understand 
localization problem in two-dimensional (2D) electron 
systems. According to the scaling theory of localiza- 
tion ||l|,^, there is no metal- insulator transition in non- 
interacting two-dimensional (2D) electron systems. One 
consequence of the theory is that there is no 2D quan- 
tum percolation transition. On the other hand, it is well 
known that there exists 2D classical percolation tran- 
sition. Thus, it is interesting to understand that in a 
realistic electron system, where electron phase coherence 
length is finite, whether the system is quantum like, or 
classical like, or even more, unlike either, belongs to a 
new class. 

In this paper, we study a 2D quantum percolation 
model [H with finite phase coherence length. The de- 
phasing mechanism is introduced by attaching current- 
conserving voltage leads to the system, a method widely 
used for one-dimensional models in mesoscopic commu- 
nity [^j^]. The conductance in such a system is derived 
and is found to contain both quantum and classical con- 
tributions. The conductance g for a finite size system is 
calculated. Through the finite size scaling analysis, we 
show clear evidence of metal-insulator transition (MIT). 
The MIT is driven by a novel type percolation transition 
which is semi-quantum and semi-classical in nature. We 
call it meta-percolation. Finally, we discuss the possible 
relevance to the newly observed 2D MIT at zero magnetic 
field. 




FIG. 1. Lattice system of our study. • denotes a lattice 
site where a phase-breaking voltage lead is attached. 

Our quantum percolation model |^ is based on a tight 
binding Hamiltonian of a 2D square lattice. 



^ = E ^-^tc^ + E ^^t^i + c+c,), (1) 
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where < ij > denotes a pair of nearest-neighbor sites and 
ei is the on-site energy. We choose et to be random, rang- 
ing from ~W/2 to W/2 to model an on-site disorder. The 
nearest-neighbor hopping matrix elements Uj is a random 
variable which assumes the values 1 (connected bonds in 
Fig.l) or (empty bonds in Fig.l) with respective proba- 
bilities P and 1 — P. In order to introduce phase-breaking 
mechanism, we attach current-conserving voltage leads 
at random lattice sites (• sites in Fig.l). The prob- 
ability where a lattice site is attached with a voltage lead 
is denoted and the hopping element between a volt- 
age lead and a lattice site is denoted ty. The sole role 
of the voltage leads is to randomize the phases of incom- 
ing and outgoing wave-functions at the lead sites while 
maintain current conservation. Using Keldysh Green's 
function formalism |^,^, it is straightforward to derive 
the multi-lead current-voltage relation: 
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where 

c^y = -J-^^ [{TjaG''o,aTai)pi{T,aG'i^Tc,j)pj] i ^ j (3) 
and 

cr„; = - ^ CTy (4) 

In the above equations, a denotes the vector space of the 
lattice system, Tai is the transmission matrix, Gq^q^ is 
the retarded (advanced) Green's function of the lattice 
system, and pi = Imgf are the density of states matrices 
for the voltage leads and measurement leads L and R, 
with being the advanced Green's function at ith lead. 
In order to satisfy current conservation, we require that 
the total current through each voltage lead to be zero, 
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namely, /j = when i ^ L, R. This restriction fixes 
the values for Vi. Under these conditions, we obtain the 
conductance between L and R 



9 = <R + ^iR- (5) 
t conductance 

L and j = R, and cr|,jj is the indirect 



Here cr^^j is the direct conductance which can be obtained 
from Eq. (3) with i 
conductance, 
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The (t];,^ is the conductance of the conductance network 

obtained from the quantum lattice model. We should 
mention that although the expression for the direct con- 
ductance ^ is identical for the pure quantum case with- 
out phase-breaking voltage leads, its value is different 
with voltage leads attached since the Green's functions 
in the formula are affected by the presence of voltage 
leads. 

We have calculated the total conductance g numeri- 
cally. In the following we present the results and discuss 
their implications. There are several parameters in our 
model as mentioned before: P (probability for Uj = 1), 
Py (probability for a voltage lead), ty (hopping between 
a voltage lead and the lattice system), W (on-site dis- 
order), and E (electron energy). The results presented 
below are for t„ = 0.1, = 2 and E — 0. We have 
done systematical studies by varying these three param- 
eters, including zero on-site disorder with W = Q, and we 
found that the qualitative results remained unchanged. 
In Fig. 2 (a) we show g as a function of P for a pTirc quan- 
tum case with P„ = 0. The calculation is done on an 
L X L square lattice. The different curves arc for differ- 
ent L's with L = 10 (solid line), L = 20 (dashed line) and 
L = 30 (long-dashed line). The important message from 
this plot is that for a given P, g decreases with increasing 
system size L, indicating an insulating behavior. This re- 
sult is consistent with the conclusion of scaling theory of 
localization, namely, in a disordered quantum system, all 
states are localized. 

In Fig. 2(b) we plot the same curves as in Fig. 2(a), ex- 
cept with Py = 0.2, i.e., 20% of lattice sites are attached 
with phase-breaking voltage leads. The striking differ- 
ence between the two plots is that in Fig. 2(b), all curves 
cross at a single point Pc — 0.62. For P > Pc, g increases 
with L, while it decreases with L for P < Pg. Thus, the 
system is metallic above Pc and insulating below Pc with 
a MIT at P = Pc. At Pc, 5 is independent of L, signify- 
ing a divergent length ^ there. Therefore, the MIT is a 
second-order phase transition. 




FIG. 2. (a) Conductance g as a function of the probability 
P for a pure quantum system with size L = 10, 20 and 30. 
(b) Same as (a), except with P„ = 0.2. Inset: Plot g vs. P at 
the vicinity of Pc for L = 10, 15, 20, 25, 30, 35 and 40. 

What are the nature of the transition and the phys- 
ical meaning of ^? The phase coherent length de- 
fines a cut-off length scale beyond which the classical 
laws for calculating a conductance are valid. One can 
divide a system into blocks with area x L^. Within 
each block, the conductance has to be determined quan- 
tum mechanically, i.e., localization physics plays an im- 
portant role. The total conductance of the system is 
obtained by applying the classical laws to the block con- 
ductances. This defines a random conductance network 
problem. The localization length Xioc of the 2D quan- 
tum system, which is always finite, determines the av- 
erage conductance and conductance fluctuations of each 
block. For the quantum percolation model studied in 
this work, A;oc depends on the probability P for a fixed 
on-site disorder W. At a small P such that Xioc *C L^, 
the block conductance guock ^ and the system 

is an insulator. With increasing P such that Xioc — L^, 
the average block conductance guock is of the order /h. 
The detail distribution of guock is determined by the cor- 
responding quantum Hamiltonian. For the quantum per- 
colation model, guock = if inside the block there is no 
at least one connected quantum path within which every 
nearest-neighbor hopping tij = 1. Hence, the distribu- 
tion function of guock always has a finite ^-function like 
peak at guock = 0, while it is a continuous function for 
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guock > 0. This defines a classical percolation problem 
in which guock may be zero or falls into a continuous 
distribution function. In this new conductance network, 
there exists a characteristic length ^, which is the aver- 
age size of the connected cluster. ^ diverges at P — P^, 
responsible for the fix point in Fig. 2(b). 
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FIG. 3. (a) Scaling plot of g vs. L/^ for L = 10, 12, 15, 
18, 20, 20, 25, 28, 30, 32, 35, 38, 40 and 45. Inset: vs. the 
probability P. (b) Plot j3 function vs. g. Inset: j3 vs. logg 
near /3 = 0. 

In Fig. 3 we show the scaling properties of our data. 
We use the standard finite-size scaling analysis The 
conductance curves for system sizes L and different P's 
can be scaled to the metallic and insulating branches. 
Those curves with P > Pc {Pc = 0.62) on the metallic 
side are scaled along the L-axis so that they coincide 
with a particular curve chosen here for P = 0.70. £,{P) is 
determined from curve shifting. The similar procedure is 
also applied to insulating side, in which case we scaled the 
curves to coincide with the curve for P = 0.54. We find 
that all the conductance g for different size L (ranging 
from 10 to 45) can be collapsed in a two-branch scaling 
function, g{L) = f{L/£,{P)), where ^{P) diverges at Pc- 



Fig. 3(a) shows the scaling function and £,{P) is shown in 
the inset. We find that ^ = iP-Pd"" with v = 1.4±0.2. 
On the insulating side, to a very good approximation. 
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Another important scaling property is the (3 function, 
defined as 



Pig) = 



d\ogg 
dlogL 



(7) 



Fig. 3(b) shows f3{g). (3 is negative for small g and 
approaches a positive value for large g. In the inset, 
P is plotted as a function of log 5 and we find that 
(3 oc \og{gJ gc) around /? ~ 0, in common with other kinds 
of MIT H. 
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FIG. 4. P„ versus Pc. 
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In our model, the phase coherence length is con- 
trolled by the probability of voltage leads Pv and the 
hopping between the a voltage lead and lattice system tv 
Although it is difficult to determine the explicit relation 
between and P^ and t^, it is evident that for a fixed t^, 
decreases with increasing P^ and approaches infinity 
(pure quantum system) at P^ = 0. We have carried out 
calculations for several values of Py for a fixed t.y = 0.1 
to determine the metal-insulator transition points. Fig. 4 
presents the resulting MIT points Pc versus Py. 

In the following, we make some remarks concerning our 
results. 

(i) We find that at the transition point Pc the critical 
conductance gc is of the order /h. In our model gc de- 
pends on the on-site disorder W. At VF = 0, — 2e^//i, 
the conductance for one conducting channel through the 
system with two spin degeneracy Thus, we associate 
the name meta-percolation to the conduction process at 
the transition point. 

(ii) The meta-percolation proposed here is different from 
the conventional classical percolation transition. Close to 
the classical percolation point, the finite-size conductance 
scales like g{L) = L~^f{L/£^) with non-zero exponent x 
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| [10{ . Thus, g{L) for different L will never cross at a 
single point as for the case of meta-percolation shown in 
Fig.2(b). 

(iii) The meta-percolation transition proposed in this 
work happens at zero temperature as long as the phase 
coherence length is finite. However, for a realistic 
electron system where dephasing is caused by electron- 
electron or electron-phonon interaction, ~ x~p which 
diverges at T = 0, and in that case, the meta-percolation 
MIT can only occur at finite temperatures. 

Before ending, we comment on the possible connection 
between the meta-percolation and the newly discovered 
MIT in 2D electron systems It is well known 

that the dielectric constant of a 2D electron system be- 
comes negative at low density where electron-electron in- 
teraction is important. A consequence of the negative 
dielectric constant, we suggest, could be the formation 
of a droplet state. The droplet state is a two-phase co- 
existence region of high density electron liquid and low 
density electron "gas". Here we call it "gas" purely for 
the reason that its density is low. In fact, in the pres- 
ence of impurities, the "gas" phase is a disordered Wigner 
crystal or a Wigner glass. The liquid phase has a better 
local conductivity compared to the "gas" phase because 
of its higher density. Thus, in the lattice model we have 
studied, the liquid phase corresponds to the region with 
nearest neighbor hopping tij = 1 while the "gas" phase 
corresponds to Uj — 0. At finite temperatures, electron 
phase coherence length is finite, thus the system is semi- 
quantum. In such a system, as we have shown, the meta- 
percolation can induce MIT and recently discovered 2D 
MIT may belong to this class. Unlike the true quantum 
phase transition, this kind of MIT disappears at T = 
since electron phase coherence length is infinity there. 
However, if dephasing mechanism can be maintained at 
T = 0, for instance, by adding magnetic impurities in the 
system, then the meta-percolation MIT can survive even 
at zero temperature. 

In summary, we have demonstrated that there is a 
novel type percolation (meta-percolation) in two dimen- 
sional systems. The meta-percolation can induce a metal- 
insulator transition and the transition is semi-quantum 
in nature. 

The work is supported by DOE under the contract 
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